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We consider properties of critical points in the interacting boson model, corre- 
sponding to flat-bottomed potentials as encountered in a second-order phase tran- 
sition between spherical and deformed 7-unstable nuclei. We show that intrinsic 
states with an effective /3-deformation reproduce the dynamics of the underlying 
non-rigid shapes. The effective deformation can be determined from the the global 
minimum of the energy surface after projection onto the appropriate symmetry. 
States of fixed N and good 0(5) symmetry projected from these intrinsic states 
provide good analytic estimates to the exact eigenstates, energies and quadrupole 
transition rates at the critical point. 

"During these moments of abstraction he seemed more intimately absolved, in 
the sense of being linked anew with the universe. 
Giuseppe de Lampedusa, "The Leopard" 

1 Introduction 

In the days that one of us (JNG) was a graduate student, group theory was con- 
sidered almost a dirty word in the nuclear physics community even though Wigner 
had introduced spin-isospin symmetry (SU(4)), Elliott had exploited the symmetry 
of the harmonic oscillator to link collective motion and the shell model (SU(3)) 
and the symmetry of the hadrons had been discovered (SU(3) again). Francesco 
Iachello changed that attitude and brought group theory front and center in nuclear 
physics and in other fields of physics. 

Franco Iachello is a descendant of a noble Sicilian family similar to that por- 
trayed in Giuseppe de Lampedusa's classic Italian novel, "The Leopard". Set in 
Sicily in 1860 at the time of the campaign for the unification of Italy, the hero, 
the Prince, struggles with how to keep the old while embracing the new. The 
Prince often took solace from the turmoil of the real world by studying astronomy 
and mathematics much the same as Franco has by his significant contributions to 
physics and group theory. 

2 Critical Points in the Geometric Collective Model 

Recently Franco has been studying two critical points associated with shape phase 
transitions in nuclei within t he g eometric framework of the collective model for 
infinite square well potentials . This model involves a Bohr Hamiltonian which 
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Table 1. Excitation energies (normalized to the energy of the first excited state) and B(E2) values 
(in units of B(E2; 2+j — > 0^ ) = 1) for the E(5) critical point^ for several N=5 calculations and 

for the experimental data of 134 Ba^3 The finite-N calculations involve the exact diagonalization 
of the critical IBM Hamiltonian (H cr i) [Eq. 1171 1. r-projected states for H cr i [Eqs. 151. 1121 with 
y = 0.314], the (7(5) limit [end] and the 0(6) limit [(A/A)(N - cr)(N + a + 4) + Bt(t + 3)]. 
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describes the dynamics of a macroscopic quadrupole shape via a differential equa- 
tion in the intrinsic quadrupole shape variables beta and gamma. In the current 
contribution we shall discuss the critical point (CP) for a second order shape phase 
transition between spherical and deformed 7-unstable nuclei, which Franco ca lled 
E(5)0 An empirical example of such a critical point has been found in 134 BaEEl 
and possibly in 104 Ru 5 , 102 Pd©and 1O8 Pd0. 

In the geometric approach the E(5) eigenfunctions ^1 are proportional to Bessel 
functions of order r + § and the corresponding eigenvalues are proportional to 
(x£ T ) 2 . Hamiltonians that are 7- unstable have an 0(5) symmetry and r is the 
0(5) quantum number. Furthermore 

X£ tT is the £-th root of these Bessel functions. A portion of an E(5)-like spectrum 
is shown in Fig. 1. It consists of states, L^ T , arranged in major families labeled by 
^ = 1,2,... and 0(5) r-multiplets (r = 0, 1, . . .) within each family. The angular 
momenta L for each r-multiplet are obtained by the usual 0(5) D 0(3) reduction 
^. The E(5) CP leads to analytic parameter-free predictions for energy ratios and 
B{E2) ratios which persist when carried over to a finite-depth potential ^ As 
seen in Table 1, the E(5) predicted values are in-between the values expected of a 
spherical vibrator [C/(5)] and a deformed 7-unstable rotor [0(6)]. 



3 Wave Function Ansatz for the 7-Unstable Limit of the 
Interacting Boson Model 

The Interacting Boson Model EH (IBM) was one of the great achievements of Akito 
Arima and Franco Iachello. This model describes low-lying quadrupole collective 
states in nuclei in terms of a system of N monopole (s) and quadrupole {d) bosons 
representing valence nucleon pairs. The IBM Hamiltonian relevant to the critical 
point of the phase transition between spherical and 7-unstable deformed nuclei has 
0(5) symmetry UJ. Its energy surface, obtained by the method of coherent states 
is 7-independent and exhibits a flat-bottomed behavior in (3 which resent- 
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Figure 1. An E(5)-like spectrum of states labeled by L^ tT . Shown are the transitions whose 
B(E2) values are given in Eq. 1121 , The E2 rates for other At = 1 transitions (not shown) are 

governed by 0(5) symmetry. Specifically, B(E2; L~^ 2 > ^li) f° r L = 4,2 are in the ratio 1 : 1 

respectively, B(E2; L+ 3 > 4+ 2 ) for L = 6,4, 3 and B(E2; L+ 3 > 2+ 2 ) for L = 4, 3, are in 

the ratios 1 : 10/21 : 2/7 : 11/21 : 5/7 : 1, respectively. Taken from refE3 



bles the infinite square-well potential used to derive the E(5) CP in the geometric 
approach Calculations with finite N values (N=5 for 134 Ba) have found that 
this critical IBM Hamiltonian can replicate numerically the E(5) CP and its ana- 
lytic predictions l ^ l ^ lZl i n t, ne present contribution, based on recent work ED we 
examine the properties and conditions that enable features of E(5) CP to occur 
in a finite system described by the interacting boson model. For that purpose we 
propose wave functions of a particular analytic form, which can simulate accurately 
the exact IBM eigenstates at the critical point. These wave functions with fixed 
iV and good 0(5) symmetry are used to derive accurate estimates for energies and 
quadrupole rates at the critical point without invoking large- N approximations. 
The proposed wave functions can be obtained by projection from intrinsic states 
with an effective /3-deformation. 

In the IBM, the 7-unstable transition region is modeled by the Hamiltonian 

H = eh d + - A A[<# -d) -^f][H.c] (1) 

with e and A positive parameters. Here ha is the d-boson number operator, H.c. 
stands for Hermitian conjugate and the dot implies a scalar product. In the Z7(5) 
limit (A = 0), the spectrum of H is harmonic, end, with = 0, 1,2 . . .N. The 
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eigenstates are classified according to the chain U(G) D U(5) D 0(5) D 0(3) 
with quantum numbers | N, rid, t, L) (for r > 6 an additional multiplicity index 
is required for complete classification). These states can be organized into sets 
characterized by rid = t + 2k. States in the lowest-energy set (k = 0) satisfy 



P \N,n d = T,T,L) =0 
Pi = d) ■ S . 



(2) 



Other sets (k > 0) are generated by | N,rid,r, L) oc (Pj ) k \ N — 2k, rid — t,t,L). 
In the 0(6) limit (e = 0), the spectrum is \A{N - a)(N + a + 4) with a = 
N, N — 2, N — 4,...0 or 1. The eigenstates are classified according to the chain 
C/(6) D 0(6) D 0(5) D 0(3) with quantum numbers | N,<j,t,L). The ground band 
has a — N and its members satisfy 



Pi | N, a = N, t, L) = 

P\ = [rft.rft _ ( s t)2] . 



(3) 



The remaining bands with a = N — 2k are generated by | N, a, r, L) oc (pj )) k \ N — 
2k,a,r,L). These results suggest that in-between the U(5) and 0(6) limits, we 
consider a ground band (£ = 1) for the Hamiltonian Q determined by the condition 



P V \Z= l;y,N, T ,L) = , 

Pt=[dt. rf t_ y(s t)2]. 

In the f/(5) basis these states are 



|e = i;y,iV,r ) i) = ^i[i + (-i)^- T ] £„ d , T |iV, n d ,r,L) , 



(4) 



(5) 



and the rid summation covers the range t < rid < N . The coefficients £n d , T have 
the explicit form 
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G^a\y) is an odd function of y, G { a'(—y) = —G)x'(y), and satisfies the follwing 
recursion relation 
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Members of the first excited band (£ = 2) have approximate wave functions of 
the form 



\£=2;y,N,T,L)=Af0PUZ = l;y,N-2,T,L) 



2{2N + y z + l) + A{y z -1)S^ T 



(JV-2) 



-1/2 



(8) 



where is defined in Eq. (|10|l below. 

The states of Eqs. (jSJ and ||BJ) have fixed N, L and good 0(5) symmetry r. 
Henceforth, for reasons to be explained below, they will be referred to as r-projected 
states. Diagonal matrix elements of the Hamiltonian Q in these states, denoted 
by E^ T = (£; y, N, t\H\^; y, N, r), can be evaluated in closed form 



N-sl 1 ? +-A(i- y ys [ 



(TV) 



E i=2 ,r = e{/V - 2Nj [2y 2 + (2N + 7y 2 - 1) + 2(y 2 1) ] } 
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If (y 2 1) C~ 2) +(V- I) 2 (2JV + tf a - 8y + 5) 



(iV-2) 



+16(i/-l)(JV + i/)5g 2) +2[(2/V + y)(2/V + y + 2) + l]} . (9) 
The E^ T are independent of L since if is an 0(5) scalar. Their expressions involve 



the quantities S^J = (£ = 1; y, jV, r|(st) fc s fc |£ = 1; y, iV, r) which are given by 
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elements of the Hamiltonian H between r-projected states in different £-bands, 
i?i,2;T = (£ = 2; y, AT, r|i?|£ = 1; y, JV, r), can be evaluated as well 

H\.2-t = 
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By techniques similar to that employed in the 0(6) limit of the IBM^, explicit ex- 
pressions of quadrupole rates can be derived for transitions between the r-projected 
states. For the relevant IBM quadrupole operator, T(E2) = S s + s^d, these tran- 
sitions are subject to the 0(5) selection rule Ar = ±1, and, as explained in the 
caption of Fig. 1, it is sufficient to focus on the B(E2) values 

B(E2; £ = 1; r + 1, L = 2r + 2 — ► £ = 1, r, L = 2r) = 
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4 Intrinsic State of the 7-Unstable Wavefunction Ansatz 



The states in Eqs. iJSJ and Eq. © can be obtained by 0(5) projection from the 
IBM intrinsic states for the ground band 



|c;iV) = (Ar!)-V2 (6 t)iV| 0; 



&t = (l+/? 2 )-l/2 



(i cosjdl + (3 sin 7 -^= (d\ + dL 2 ^ + 



(13) 



V2 

and for the (3 band respectively 

\(3;N) =MpPl\c;N-2) , (14) 

provided y — 1 . The expressions in Eqs. (5)-(12) depend on the so far unspecified 
parameter y. Normally, the equilibrium value of /?, and hence y, is chosen as the 
global minimum of the intrinsic energy surface 

determined from the expectation value of H in the intrinsic state (|13|) . This is a 
standard procedure for a Hamiltonian describing nuclei with rigid shapes, for which 
the global minimum is deep and well-localized. Such is the case for the Hamiltonian 
of Eq. whose intrinsic energy surface has the form 

E(f3) = E Q + N(N - 1) /3 2 (1 + [3 2 )- 2 [a + cf3 2 ] 

a = e-A , c = e , e = e/(iV-l) (15) 

with Eo = jAN(N — 1) a constant. The topology of the energy surface is such that 

a > global spherical minimum at j3 = 

a < global deformed minimum at (3 — a /(a — 2c) . 

When a / the intrinsic energy surface behaves quadratically (~ 1 ) near the 
single minimum, and the standard procedure of determining the equilibrium value 
of (3 from the global minimum of E{(3) is applicable. However, near the critical 
point of the phase transition an alternative procedure is required. 

5 7-unstable Wave function Ansatz at the Critical Point 

The IBM Hamiltonian, H cr i, at the critical point of the U(5) — 0(6) phase transition 
corresponds to a special choice of parameters in the Hamiltonian of Eq. Q 

H cri : e=(N-l)A, (17) 

for which a = in Eq. I|15|l . and the corresponding energy surface reduces to 

E(P) = E a + AN(N-l)p 4 (l+(3 2 )- 2 . (18) 

In this case, the intrinsic energy surface E(f3), shown in Fig. (2a), has a flat behavior 
(~ /3 4 ) for small /?, an inflection point at f3 = 1 and approaches a constant for 
large (3. The global minimum at /3 = is not well-localized and E((3) exhibits 
considerable instability in (3, resembling a square- well potential for < (3 < 1. 
Under such circumstances fluctuations in /3 are large and play a significant role 
in the dynamics. Some of their effect can be taken into account by introducing 
into the intrinsic states of Eqs. i|13|) and (|14l) an effective /^-deformation. The 
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Figure 2. Energy surfaces of the critical IBM Hamiltonian H cr i 1171 with N = 5 and A = 1. 
(a) Intrinsic energy surface E(/3), Eq. 1181 [solid line], and its approximation by a square-well 
potential [dashed line], (b) 0(5) proje cted energy surface _E^ =1 t _q (y), Eq. JHJ. The global 
minimum is at y = 0.314. Taken from refE3 



effective deformation is expected to be in the range < y — 1 < 1, in-between 
the respective U(5) and 0(6) value of (3. This will enable a reproduction of E(5) 
characteristic signatures which are in-between these limits (see Table I) . In contrast 
to E((3) of Eq. H18[). we see from Fig. (2b) that the 0(5) projected energy surface 
E^ =lT=0 (y) of H cri (Eqs. 10 and (|17|) with N ~ 5), does have a stable minimum 
at a certain value of y, which we interpret as an effective /3-dcformation. This 
procedure, based on variation after projection, is in the spirit of in which it is 
shown that in finite boson systems, a 7-unstable 0(6) state can be generated from 
a rigid triaxial intrinsic state with an effective 7-deformation of 30°. In the present 
case the 7-instability is treated exactly by means of 0(5) symmetry, while the /3- 
instability is treated by means of an effective deformation. The appropriate value of 
y can be used to evaluate the band-mixing, rj T (y) = g~^rgj ■ A small value of rj T 
will ensure that the r-projected states of Eqs. (JSJ), |(HJ form a good representation 
of the actual eigenstates of H cri , and turn the expressions of Eqs. ©, (fTSfl into 
meaningful estimates for energies and quadrupole transition rates at the critical 
point. 
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Table 2. U(5) decomposition (in %) of the T states for N = 5. The calculated values are 
obtained from the r-projected states, Eqs. JSJ, JSJ with y = 0.314. The exact values are obtained 
from numerical diagonalization of the critical IBM Hamiltonian H cr i, Eq. 1171 . 
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6 Comparison with E(5) and Experiment 

To test the suggested procedure we compare in Table 2 the C/(5) decomposition of 
exact eigenstates obtained from numerical diagonalization of H cr i for TV = 5 with 
that calculated from the r-projected states with y = 0.314 [the global minimum 
of Eifl(y)}. As can be seen, the latter provide a good approximation to the exact 
eigenstates (the corresponding band-mixing is r\ T = 0.12, 3.53, 4.14, 3.05% for r = 
0, 1, 2, 3). This agreement in the structure of wave functions is translated also into 
an agreement in energies and B(E2) values as shown in Table 1. The results of 
Table 1 and 2 clearly demonstrate the ability of the suggested procedure to provide 
analytic and accurate estimates to the exact fmite-N calculations of the critical 
IBM Hamiltonian, which in-turn agree with the experimental data in 134 Ba and 
captures the essential features of the E(5) critical point. 



7 Large N Limit 

In the large N limit, using Stirling's Formula in Eq. ©, we obtain 

G'"b)^(^ 2 "^ (2p)!!(2 ;+ 2 „ +1) ,, . 



where x = Ny. Therefore, 



si N T \y)^s k>T (x) , (20) 



that is, Sl N J(y) becomes a function of x only. This suggests that the energy spec- 
trum depends only on x and not N and y separately as N goes to infinity, except 
for possibly an overall scale. To test this in Fig. 3 we plot the energy ratio at the 
critical point 

R= E^ T=2 -E^, T=0 (21) 

£j£ = l,T=l — £<f=l,T=0 

for several values of N as a function of x. Indeed as N increases, the curves asymp- 
tote to one universal curve. This means that the quadrupole beta deformation of 
the geometrical model is proportional to y/~N (3. This observation may be useful in 
relating the predictions of E(5) CP to the IBM. 
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1 2 3 4 5 

x 

Figure 3. The energy ratio R at the critical point defined in Eq. 1211 as a function of x = Ny, 
for different values of N. Note that the vertical axis is displaced from zero. 



8 Summary and Future Outlook 

In thi s co ntribution we have examined properties of a critical point in a finite 
system^! ^e have focused on the E(5) critical point relevant to a second-order 
shape phase transition between spherical and deformed 7-unstable nuclei. At the 
critical point of such a phase transition the intrinsic energy surface is flat and there 
is no stable minimum value of the deformation. However, for a finite system, we 
have shown that there is an effective deformation which can describe the dynamics 
at the critical point. The effective deformation is determined by minimizing the 
energy surface after projection onto the appropriate symmetries. States of finite 
N and good 0(5) symmetry, projected from intrinsic states with this effective 
deformation simulate accurately the exact eigenstates, and can be used to derive 
analytic estimates for energies and quadrupole transition rates at the critical point. 

In the future we shall explore in depth the ./V-dependence as well as the large- N 
limit and its relationship to the geometric E(5) critical point. We shall also explore 
the first-order critical point in the phase transition from a spherical vibrator and 
an axially symmetric deformed nucleus within the IBM. 
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